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Abstract 

. . . The search for a quantum theory of gravity has become one of the most well-known problems in 
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^^ ■ led to a search for a new theory of gravity that, while still agreeing with measured observations, is 

^r< renormalizable. In this paper, I show that, given a "vortex" model of elementary particles in which 



< 
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rest mass derives from intrinsic spin and polarization, a Yang-Mills force with a S0(4,l) group sym- 
metry predicts post-Newtonian N-body motion such as solar system observations of gravitational 



i-Q ■ behavior as well as binary pulsar precession and orbital speed-up caused by gravitational radiation- 
df 

r~| ■ reaction. I give a definition of the Yang-Mills theory on a lattice graph such that a background 
-)— > ■ 

■ ■ manifold is not required to define the theory. Using a homogeneous, isotropic universe model, I 

show that this theory does not contradict cosmological observations of Type la supernovae, Baryon 

►-j^ ■ Acoustic Oscillation, and the Cosmic Microwave Background. Most importantly, it agrees with 

QQ i the accelerating expansion of the universe as a consequence of the de Sitter group Lie algebra — an 

Tij- ■ acceleration that I show does not occur in the Poincare group approximation — suggesting that 



the de Sitter group symmetry explains dark energy. In addition, because it is a generic massless, 
semi-simple Yang-Mills theory, it is mathematically proved to be a perturbatively renormalizable 
quantum theory. 
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I. INTRODUCTION 

That general relativity is able to subsume Newtonian gravity and explain phenomena 
that do not fit into the Newtonian framework such as Mercury's perihelion precession, light 
deflection, and gravitational red-shift as well as being compatible with special relativity has 
brought it wide acceptance. In recent years, however, as measurement tools have grown more 
accurate and new observations made, the necessity of introducing either tunable parameters 
or exotic forms of matter and energy to make general relativity fit those measurements has 
brought the theory into question. Observations of galactic rotational curves and gravitational 
lensing [1] [2] [3] have demanded the introduction of dark matter and accelerating expansion 
of the universe [4] dark energy. At present there is no consensus on what these substances 
are, whether they even exist, or whether they are features of a more complete theory of 
gravity. 

Besides the problems with general relativity at the macroscopic scale, achieving a quan- 
tum theory of gravity has become one of the most significant unsolved problems in physics. 
Attempts to place the Einstein-Hilbert action, 

Seh = / d'^Xy/^R, (1) 

where R is the Ricci scalar, into a functional path integral ensemble, have all failed, giving 
nonsense results. The immediate source of the difficulty is that the theory is not renormaliz- 
able creating an infinite number of counterterms in perturbation expansions. Another, less 
theoretically troubling but mathematically difficult issue is that the expansion of the La- 
grangian fails to terminate because of the volume element ^J~^ and the inverse metric g^" . 
Therefore, unlike the actions of other forces, the gravitational action is not finite polynomial. 
A similar disaster occurs quantizing the weak force but disappears when it is unified with 
the electromagnetic force, motivating the quest to unify gravity with the other forces in the 
hope that it will become renormalizable [5]. This unified theory, however, would not solve 
all of gravity's problems because, by design, it would become general relativity at galactic 
and cosmological scales where dark matter and energy have an effect. Therefore, a unified 
theory would not necessarily solve these large scale, classical problems without additional, 
low-energy features. 

Although elegant and surprisingly accurate until recent decades, because the curved 
spacetime assumption causes such deep, unresolvable problems at the quantum level, there 



is some motivation to find an alternative to spacetime curvature as the source of gravity. 
Despite its current acceptance, curved spacetime was a counterintuitive choice for a theory 
of gravity given the success of electromagnetism in the 19th century. As Wald points out 

[6], 

Maxwell's theory is a remarkably successful theory of electricity, magnetism, and 
light which is beautifully incorporated into the framework of special relativity. 
Therefore, one might expect that the next logical step would have been to develop 
a new theory of the other classical force, gravitation, which would generalize 
Newton's theory and make it compatible with special relativity in the same way 
that Maxwell's theory generalized Coulomb's electrostatics. However, Einstein 
chose an entirely different path[.] 

Because curved spacetime has never been precisely measured, the path not taken still lies 
open. 

Indeed, we can ask how precisely general relativity has been confirmed, and, as this 
paper explores in detail, the answer is two-fold: (1) if general relativity's field equations 
were as simple as Maxwell's equations or not very much more complicated, then it could 
be considered a well-confirmed theory because its linearized version has a great deal of 
evidence behind it, and (2) although observation has confirmed many effects, many of the 
highly non-linear portions of the theory have been unobservable in any precise way because 
of the weakness of gravity. All precision experiments (classic tests) of general relativity 
have been done within the solar system (with the exception of binary pulsar precession) 
where spherically symmetric, weak gravity prevails and only the first order post-Newtonian 
Einstein-Infeld-Hoffman (EIH) equations have been confirmed [7] [6]. An example of the best 
recent evidence for strong field general relativity is the periastron precession of the double 
pulsar system PSR J0737-3039A/B [8] , which is a first post-Newtonian order confirmation. 
Furthermore, recent measurements of the cosmic microwave background, where spacetime 
curvature is most likely to appear, show none at all [9] with the infiationary theory as a 
potential explanation [10]. 

The difference between the EIH equations and the full equations of GR has invited a 
large number of alternative theories of gravity such as are discussed in [11], but virtually all 
are metric theories and retain the curved spacetime picture with the associated problems 



mentioned above. Non-metric theories, by contrast, offer the hope of ehminating gravity's 
problems with quantization at the cost of also eliminating the elegant curved spacetime 
approach. It has been suggested that the Standard Model approach, where forces have 
Yang-Mills actions, may be more suitable for a quantum theory of gravity because of its 
success at explaining the other forces. The tetrad or veirbein formulation of the Einstein- 
Hilbert action, for example, has a structure similar to the Yang-Mills [5]. Conformal gravity 
has a Yang-Mills structure as well but tends to diverge strongly from Newtonian gravity 
over large length scales [12]. None of these theories have achieved significant success in 
making predictions that general relativity is not capable of making nor in resolving, entirely, 
the problems with the quantum theory. Thus, there has been, as yet, no motivation for 
replacing general relativity with any of these nor of accepting any as quantum gravity's true 
representation. 

In this paper, I present a Yang-Mills theory of gravity with a S0(4,l) or de Sitter group 
without a background metric. Because a non-Abelian S0(4,l) theory has local Lorentz sym- 
metry, the theory complies with the Einstein Equivalence Principle (EEP) and the Strong 
Equivalence Principle (SEP) (which even most alternate metric theories do not satisfy). 
All physical theories are Poincare symmetric, but, given that the Poincare group is not 
semi-simple, the de Sitter group is a natural choice of a semi-simple extension. This paper 
represents the first time that a Yang-Mills theory based on this group has been shown to 
agree with precise gravitational observations. 

Yang-Mills theories in general have a finite, polynomial action with a renormalizahle 
quantization provided the coupling constant has non-negative mass dimension. Since G has 
a mass-dimension of —2, we turn to an unknown non-dimensional gravitational constant, a^, 
for quantum predictions. Because the coupling constant is non-dimensional, mass, energy, 
momentum, and angular momentum are all quantized, which has already been suggested in 
the context of other theories of gravity [13]. Indeed, any quantum field can be decomposed 
into finitely many modes, 

M/2 

<P= Y. «ne*'""^ (2) 

n=-M/2 

which constrains the number of degrees of freedom of the field to be a finite number, M. 
The quantity M is simply a regularization (energy cutoff) which is taken care of by renor- 



malization, and, if the kn are closely spaced then, we can take the continuum approximation, 

I'M/2 ,,^ ^ 

<Pi^,) ^ / j^TYAkne''''^^. (3) 

J-M/2 {^'^r 

This is an approximation only however, similar to the thermodynamic limit approximation 
of classical fluids, taking the number of molecules to infinity, and, at a fundamental level, 
the number of degrees of freedom of the field remains finite. 

With a lattice method, I show that the background coordinate system normally present 
in Yang-Mills theory can be replaced with an artificial, "bookkeeping" coordinate system 
where all distances, time intervals, and torsions in spacetime are derived from the Yang-Mills 
potential. In this way, the Yang-Mills action is the same as the standard Yang-Mills action 
for SU(N) theories, but its potential alone is the metric of spacetime. 

Rest mass is derived from the energy of spin and polarization and this is the second 
component to the theory. In any frame, these vectors form an anti-symmetric tensor S^'^ 
such that polarization is p* = 5*°* and spin is s* = e^i^S^^ where summation is implied. Spin 
and polarization have a conserved current (similar to the symmetric electromagnetic stress 
energy tensor [14]), 



rpfMl/ 



S^'^Sx" + -v^^'SapS''^ 



(4) 



where r]'^'^ = diag(— 1, 1, 1, 1) is the Minkowski metric. This tensor is both symmetric and 
traceless. This implies that mass derives relativistically from an intrinsic kinetic energy of 
spin and oscillation. 

With these definitions, this Yang-Mills theory, called vortex gravity because of the vortex 
particle interpretation that spin and oscillation implies, predicts N-body observations such 
as pertain to the Solar System and binary pulsars including radiation reaction, all with a 
finite polynomial, positive definite action that is renormalizable. 

I show that the theory offers a cosmology slightly different from the ACDM dual epoch 
model. In this cosmology, I make a distinction between Doppler and gravitational redshift, 
i.e. redshift caused by matter moving within space and redshift caused by the expansion 
of space itself. For the expansion of matter, I derive a linear coasting cosmology. The 
linear cosmology has been shown to fit observations such as Type-la supernovae and the 
Cosmic Microwave Background related to the later universe well [15] [16]. In additional, I 
derive a gravitational potential which determines the degree of gravitational redshift. This 
leads to an interesting result: while in the Poincare group approximation of the de Sitter 



theory, the hnear coasting model dominates the later universe with the gravitational redshift 
vanishing, in the de Sitter theory the combination of the Doppler and gravitational redshift 
predicts a linear acceleration in total observed redshift. I explain this phenomenon by 
observing that, in the de Sitter Lie algebra, momentum produces spin by the coupling of 
momentum generators [V, V] = iAf^" where V are momentum generators and Af^" are 
spin generators. In the Poincare group spin can only be produced if spin is already present 
because [V'^, V] = 0. The ability of the de Sitter theory to produce spin from momentum 
without a source of spin causes the expansion of the universe to produce gravitons that 
couple to spin which in turn create gravitons that couple to momentum. The back and 
forth production of gravitons leads to an accelerating redshift. This explanation is similar 
to the explanation for Thomas precession where the orbital motion of an electron increases 
its spin — an anomaly explained in the early 20th century which can be shown from the 
Lorentz Lie algebra [14]. Hence, "dark energy" can be construed to be gravitons produced 
because of the relationship between momentum and spin in a de Sitter symmetry group. No 
cosmological constant is required, and no vacuum energy is modelled. The effect is produced 
directly from the Yang-Mills equations. Observations currently attributed to dark matter 
such as galactic rotational curves and anomalous lensing are not addressed in this paper. 

The paper is organized as follows: Sec. II A describes the conserved quantities of mass 
and linear momentum, oscillation momentum, and angular momentum, how these are rep- 
resented in the theory, and how they appear intrinsically through the vortex model; Sec. 
II B derives the Yang-Mills theory with only a lattice and the gravitational potential and 
no background manifold. This section also explains how the bookkeeping coordinate system 
works; Sec. II C shows a derivation of the field equations for the S0(4,l) Yang-Mills theory 
and its formulation as two covariant field equations; Sec. HE derives the equation for parti- 
cle motion and shows that it is the same as the geodesic equation; Sec. II F gives a discussion 
for how the theory satisfies the Strong Equivalence Principle (SEP) which GR also satisfies; 
Sec. V discusses novel predictions made by the theory; Sec. II G briefly mentions the mathe- 
matical justification for a renormalizable quantization; Sec. Ill shows that the theory agrees 
with measurements within the solar system including redshift, perihelion precession, time 
dilation and light bending (IIIB), binary pulsar precession (IIIC), the quadrupole formula 
for gravitational radiation as well as showing that there are no dipole/monopole moments 
(III CI), and the equations for an expanding Robertson- Walker universe (HID); Sec. IV 
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explores any possible experimental disagreement with the theory; Sec. V discusses tests that 
may confirm or deny the theory's predictions; and Sec. VI discusses related work. There 
are also two Appendices. 

II. THE THEORY 

Although ad hoc gravitational theories exist, any strong theory tends to be based on 
a guiding principle. In developing this theory, the guiding principle is that gravity is an 
ordinary force in the standard, SU(3)xSU(2)xU(l) model. This does not imply, as is 
frequently mentioned in the literature, that gravity needs a background manifold to be 
defined as a non-metric theory. I show that it does not provided one accept a lattice graph 
construction to spacetime. In some sense, the lattice graph is no more presumptive in a 
discrete sense than Einstein's differential manifold and point-coincidence assumptions are in 
a continuous sense [17]. They both simply assume that events take place at points and that 
points have neighboring points. Indeed, the lattice version is the weaker of the two because 
it does not require continuity meaning that neighboring points need not be measurably 
close. An essential feature of the following theory is that all the physical requirements are 
supported in the literature, and all have physical explanations based on well-established 
principles derived from experiment. 

A. Conserved quantities of the vortex model in the classical limit 

The vortex model of elementary particles was first introduced by Lord Kelvin for atoms 
and since expanded with modern intrinsic quantum properties of particles [18]. Even pre- 
dating the development of special relativity, the mass of the electron has been presumed to 
derive from its electrical energy [19], and hence the stress-energy tensor of an electron, and, 
by extension to quarks and other particles, all matter, must, classically, take the general 
form of the Maxwell stress-energy tensor. For neutral particles, their stress energy is based 
on spin and oscillation rather than magnetic and electric fields, but the critical point is that 
spin/oscillation has, not only mass, but kinetic energy in the form of static pressure. In 
other words, particles cannot be represented as "point masses" , blobs of matter where the 
only non-zero component of their stress-energy tensors is T^^ = m (in the rest frame). The 



stress-energy tensors have space-space components as well, 5jkT^^ = m. Because they are 
positive, these space-space components cannot be cancelled out in scaling up to macroscopic 
bodies. 

The vortex model proposes that all mass derives from intrinsic quantum properties of 
spin and oscillation. Spin and oscillation are not separate entities, and each particle has a 
covariant, anti-symmetric spin tensor, s^" = —s'^^, forming a single 4-dimensional spin in 
the six planes of a Minkowski space. In any frame, oscillation occurs in the direction of the 
polarization vector P* = s*^* and spin in the plane with spin vector S** = e^^^s^^. 

The 5x5 matrix current of the S0(4,l) Yang-Mills theory, J^, describes the motion of 
particles (see A for an overview of the de Sitter group). This current, written in component 
form, is. 



J^. = T.^iy) + S^^pM'^^ 



(5) 



where S'^q,/3 represents Lorentz boosts (oscillations) and rotations (total angular momentum), 
and T^iy represents the stress energy momentum tensor which is conserved by translations, 
where /i, z/ = 0, 1, 2, 3. 

I define the conserved currents in terms of spin and polarization[35] vectors. Classical 
particles oscillate and spin at the speed of light, c = 1; therefore, we have rest frame 
polarization vector p* and rest frame spin vector s*. 

From the spin/polarization vectors, s* and p*, a tensor is formed which can then be in 
any frame. 
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s. 
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pi p2 p3 y^ 



—p^ s^ s^ 

_p2 _g3 Q gl 
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(6) 



The first current, total energy, is given in any frame by. 
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(7) 



Like the electromagnetic stress energy tensor, this has the form (in Cartesian coordinates 



{x,y,z)) [14], 
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(8) 



where U = p x s, Pij = piPj + Si-Sj - ^{jf + s^)5ij, p^ = {p^y + {pyf + ipz) % and s^ = 

{sx? + {syf + {s.y. 

Similar to the electromagnetic Poynting vector, S'poynting = \E x B, the vector U/m is a 
particle's flux. The energy flux is U. Particles in the rest frame have no energy flux. This 
gives a deflnition of a massive particle's rest frame: it is the frame where the polarization 
and spin vectors are parallel, p x s = 0. This means that in every particle's rest frame, 
the plane of oscillation is orthogonal to the plane of spin. Classically, one can visualize this 
as a point particle pursuing a periodic helical motion revolving about its axis of spin while 
oscillating normal to the plane of spin like a tiny vortex. As the particle increases velocity, 
the angle between the planes decreases until they become parallel at the speed of light. 
Therefore, the two 3- vectors are plus or minus the same direction, p* = ±s*. 

The conserved current for massive particles in the rest frame has a convenient represen- 
tation in spherical coordinates. 



Too = Trr = ip^ + s^ 
Because p^ + s^ = m (with c = 1) in the rest frame. 



Ton — i r, 



m 



(9) 



(10) 



Matter and antimatter both have positive relativistic mass (as do photons) as a natural 
consequence of the positiveness of the polarization and spin energy. This energy does not 
change its sign. For example, in a conversion experiment from a electron-positron pair in a 
head-on collision to a pair of photons emerging in opposite directions, e~ +e+ = 27, because 
the total momenta sum to zero (neglecting spins), the relativistic masses going in must sum 
to those coming out. Hence, 

n- + n; = 2T,i. (11) 

Because electrons and positrons have the same relativistic mass if travelling at the same 



speed, we have Tqq = TZ. This implies that matter, antimatter, and photons are all 



gravitationally attracted to one another and, further, that antimatter behaves the same way 
in a gravitational field as matter. (This result is already predicted by general relativity.) 

Finally, the conserved currents described here should not be confused with the definition 
of stress-energy-momentum in general relativity. General relativity provides a definition 
of energy based on the symmetries of spacetime geometry. If we have a Lagrangian that 
depends on the metric field, C{g^i,), then energy is simply the variation of that Lagrangian 
with respect to the field, 

1 ^^_ ^ j.^u_ (^2) 



^/-g dg^u 
This is not the tensor, T'^'^, we have been discussing above. General relativity assumes 

diffeomorphism-covariance which implies that the quantity 12 is conserved by Noether's the- 
orem. Kretschmann pointed out to Einstein himself (and Einstein conceded), however, that 
general covariance, i.e. coordinate covariance, was a fictitious symmetry group generating 
no physical predictions [17]. Many theories can be made that are generally covariant and 
predict different physical outcomes. The physical restrictions on the theory are important 
and fix the theory. Diffeomorphism-covariance in general relativity is an active symmetry 
that is unrestricted. This is not so in the Yang-Mills theory because the gauge prefers a 
particular coordinate system. In other words, once the gauge is fixed it specifies a spe- 
cific set of reference frames at each spacetime location. With a preferred set of reference 
frames, although the coordinate system can be changed at will, the theory is not actively 
diffeomorphism-covariant. Hence, Noether's theorem does not apply to the metric, and the 
stress-energy-momentum tensor of General Relativity is not a conserved quantity. This is 
not the same as having a preferred coordinate system overall, only that fixing the gauge 
of the Yang-Mills theory destroys the diffeomorphic gauge freedom that general relativity 
assumes. In the next section I will show this mathematically. 

B. Derivation of a Yang-Mills theory without a background metric 

Yang-Mills theory is typically defined with a Minkowski background metric, 1]^^. In this 
section, I derive a Yang-Mills theory without a background manifold using a dual labelling 
approach where I associate each vertex on a four dimensional lattice with two labels: (1)4- 
vector, i/a = {yo, yi,y2, ys), labels derived from the gravitational potential that represent real 
distances and torsions between vertices and (2) a second set of labels, Xa = {xo,Xi,X2,X3), 
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that are equidistant between neighboring vertices. The mapping Xa — )■ Va is one-to-one and 
onto since there is exactly one of each label for each vertex. At no time does a background 
need to be invoked, but, when I take the continuum limit, a metric appears. In this section 
only, I will use the beginning letters of the latin alphabet a,b,c, . . . ,hto represent 4- vectors. 
(I prefer not to use the usual Greek letters fi, v because they are typically used with a metric 
and a diffeomorphic manifold. Here we assume no general covariant structure because there 
is no manifold.) These letters range over 0,1,2,3. I will use 1/4 and 3:4 to refer to the 
additional, 5th dimension that completes de Sitter group vectors. I will use the middle latin 
letters i, j, /c, / (in this section only) to refer to entire vectors. When I do so, the vectors will 
be bolded, e.g., yi. 

Consider a four dimensional lattice and label each vertex with four numbers yo? Vxi 1/2? l/s- 
We do not yet know what these labels are. Let time be imaginary t — )■ ir. We will analytically 
continue our result back to 5*0(4, 1) later. Consider three vertices in a row |/,z, and w, let 
Uyw{z) be a small rotation matrix with e > the small parameter. Uyw = exp{eAa) can be 
split into two parts: Gab and Hate by the following: an SO (5) potential can be written as a 
sum of generators: 

Aa = GabVb + HabcMhc, (13) 

where Vb are generators for rotations in the yb-yA plane and Mbc are generators for SO (4) 
for rotations in the yb-yc plane where a,b,c = 0, 1, 2, 3. We will refer to Gab and Habc as the 
gravitational potentials. 
To linear order in e, 

Uy^=I + eAa = 1 + e{GabVb + HabcMbc). 

Each of the generator matrices is a 5 by 5 matrix. Let y = {yo,yi,y2,y3,y4) and w = 
{wq, wi,W2,W3,W4), then we define 

w := [J + e{GabVb + HabcM,,)]y. (14) 

We need to define only one point in the entire lattice, y = (0,0,0,0,1/4), in order to define 
all the other points using this procedure. 

Getting rid of the generator matrices, the equation 14 implies the relationship between 
the 4-vectors, 

Wa = eHabcyc + (^Gab + yb, (15) 
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again, all to linear order in e. The first term in the right hand side of 15 represents the 
torsion or twisting of a vector as it travels from y through z to end at w. The second term 
on the right hand side represents the translation. (Note: if Gabi"^) = then there is no 
translation and the two lattice points are effectively in the same location. This means that 
it is possible for spacetime to collapse in on itself as in the Big Bang.) 

Thus, the labels of the lattice, y, which previously had no meaning, are now defined in 
terms of the gravitational potentials by hopping from point to point outward from the only 
pre-defined point y (which is arbitrary of course). A similar procedure is used to define the 
natural numbers in terms of a successor function S, where is the only predefined point, 
from Peano's axioms. The gravitational potentials, G and H, are successor functions that 
define the spacetime lattice. 

Having defined the labels as coordinates in a discrete spacetime, we can derive Yang- Mills 
theory using Wilson's method for lattices. First, define the second labelling system by the 
same method but using trivial potentials G^^ = 6ab, the Kronecker delta such that 6aa = 1 
and 6ab = for a ^ b, and H'^^^^ = 0. This labelling system is denoted by Xa and corresponds 
to a lattice where each vertex is a distance e/2 from each of its six neighbors. 

Consider a plaquette p between four vertices in the a — b plane between Xi, Xj, x^, Xi, and 
back to Xi, and assume that each vertex is two lattice edges from its two neighbors. The 
action on the plaquette is given by Sp = Re tr UijUjkUkiUu. Now, Uij = exp(eAa(x)) where 
X = (xi + Xi)/2 is the midpoint between Xi and Xj and likewise for Ujk, Um, and Uu [5]. 



'klUli 

Now, some matrix algebra needs to be done and can be found in any standard text book 



f/.,f/,. = e-'^'-e~'^'b = e"^«+^6)+l^'[^'"^^l+o(^') 



on quantum field theory such as [5] or in Wilson's original paper [? ] and we have, 
which simplifies to 



U^,UikUklUl^ = I+e^ 

e e 



[Aa — Ag) \Ab — A^ 



b) 



+ [Aa, A 



(^ - K) {Ab - A- 



[Aa, A 



to the required order. Summing over all plaquettes and we have a constant term from the 
first which we call the vacuum energy (note that here the gravitational potential does not 
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couple to vacuum energy), the middle term cancels out, and we are left with an equation 
for the total action, 

S.^Retr (4(^1^ -(:^ + l^.. All)' 

p \ L J / 

Taking e — > gives us the Yang-Mills potential and continuing to real time from imaginary 
time changes the 6ab operators into rjat operators, 

S = T^ I (^'^xrjacVbdFabFcd, (16) 

where g is the coupling constant and 

Fat = daAb-dbAa+[Aa,Ab]. (17) 

The Minkowski metric has appeared because we continued from an SO (5) group to a 
S0(4,l) group and our second labelling system is generated with G^^ = 6ab- In real time, 
that generator becomes G^^ = rjab- 

The so-called "background" metric of our theory, then, is the metric over an artificial 
labelling system Xa, while the true structure of spacetime is described by the first labelling 
system ya- The artificial coordinate system is anything we choose it to be. It need not be 
Cartesian. The only requirement is that we have the ability to map every coordinate point 
Xa to every true coordinate point ya in a one-to-one, onto fashion. Therefore, we have the 
freedom to choose coordinate systems that are useful. 

It is important to take away from this section that the Xa coordinates, metric, and diffeo- 
morphic symmetry on the coordinates are not real in the same way that the ya coordinates, 
gravitational potentials, and locally de Sitter symmetry are real. The coordinate system, Xa, 
is generally covariant in the continuum limit, but the diffeomorphic symmetry is vacuous 
because the coordinate points are not measurable. Only the first coordinate system, ya, 
is measurable. Distances and time can be measured with rods and clocks. Torsion can be 
measured by detecting particle spin and oscillation. Hence, the coordinates ya are something 
an experimentalist can discover, but Xa can only be invented. Therefore, the S0(4,l) action 
16 describes reality while the diffeomorphic coordinate system, Xa, describes a bookkeeping 
system imposed on reality. 
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C. Yang-Mills theory and derivation of the field equations 

In this section, I begin with standard Yang-Mills theory and derive a locally de Sitter 
theory of quantum gravity using the de Sitter Lie algebra. Yang-Mills theory describes 
forces as exchanges of gauge bosons based on a group symmetry. The Standard Model of 
quantum field theory currently contains three forces, electromagnetism, the weak force, and 
the strong force in a U(1)(8>SU(2)(8>SU(3) group symmetry. 

A conventional SU(N) Yang-Mills theory has action, 

S = -^^Jd'xF,^F'^\ (18) 

where g is the coupling constant and 

F^u = d^A, - d,A^ - t[A^, A,], (19) 

is the "force" where A^ is a matrix potential from the SU(N) group, written with indexes 
as A^ab- For the S0(4,l) group these are 5x5 matrices. (Note that in the previous section 
we absorbed a factor of i into A^, to make the matrices real. The two expressions for F^^^ 
are equivalent.) The theory is invariant under a gauge transformation U{x) such that A' = 
U{x)A^W{x)—i/gd^U{x)U'^{x). For infinitesimal transformations U{x) = e*^*-^^ ~ I + ix{x) 
we have A'^ = A^, + df,x{x) + i[x{x), A^]. 

Yang-Mills equations are often expressed in group component form, where the action is, 

F;, = d,At - 9,A^ + r'^A'^Al, (21) 

and f"-^'^ is the group symbol. Each component of the potential represents a field of spin-1 
gauge bosons. 

Rather than decomposing into pure component form which is not covariant, it is more 
useful to put it into a covariant form. The SO (4,1) potential can be written as a sum of 
generators of Lorentz rotations and boosts and of boost/rotations with respect to X5, 

A, = G,,^ + H^,^M-\ (22) 

where the potential G^j^y couples to T^'^ and H^^x to S^'^^. In this sense, the theory resembles 
Einstein-Cartan theory. 
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Following the de Sitter algebra of A, the "forces" of gravity decompose into a covariant 
form with two equations analogous to the Einstein field equations and the Cartan torsion 
equations respectively, 

E^ux = d^Gi^x — dyG^x + 
V'^'^ {GfjtxHucrp — G^pH^^x — G^xHfjiap + G^pH^^^x) , (23) 

and 

dpHyafi — dyHp_al3 + Gp^aGyfj — Gp^pG^a + ^fiiya/S — ^p.u(5a (24) 



puaP 'JpJ^uaP 



where 



^pual3 — f]'^'^ {HpaaHupP — H^„aHu/3p — 

HpaaHypP + HpaaHufip) ■ (25) 

Thus, the action of gravity is, 

Sgravity = " T / d X E^^xE^'^ + Fp.yapF^'"^ . (26) 

Here the coupling constant is ag = g^. 

The classical equations of motion of any Yang-Mills theory can be found by the variation 
of the action, 6S = 0, via the Euler-Lagrange equations, 

^' V9(9M7y j "^47^°' ^^^^ 

where C = F'^^^FJ^" is the Lagrangian. Evaluating the Euler-Lagrange equations for the 
Yang-Mills Lagrangian, the equations of motion are, 

d^F^, - t[Fp,, A^] = -Sng'j,, (28) 

where J^ is the conserved current. 

Applying the definition of the group symbol again and the previous definitions 23 and 24 
to the field equations (28) yields the following equations: 

d'^E^.x + V''' {E^,xH\p - Ep,pH\x- 

G^xFpu^p + G%Fp,,x) = SttT^x (29) 

d^FpyaP + EpyaG^ ji — EpypG^a + ^uap — ^u/Sa = ^T^S^ajS (30) 
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where 

Fp,yacrH^pP + Fp^yaaH^pp) (31) 

and the usual negative sign on the couphng constant has been absorbed into the current 
in the same way that —i is absorbed into the potential. The coupling constant, Qg = g, is 
non-dimensional, but, because the size of the mass quantum is not known, I use the standard 
coupling constant, k = G/c^, in units such that k, = 1. 
Boundary conditions are: 

Gp^ -^ r]p^, dxGp^ -^ (32) 

and 

Hpux -^ 0, daHp^x -^ (33) 

as Xp — )■ oo. 

Gauge transformations can be done in component form. Consider the 5*0(4, 1) de Sitter 
transformation A{x) = ^p{x)V'^ + XfiuM'^". The gauge transformation is 

G'pu = Gpu + dpCu{x) + r]P^{^x{x)Hppy - ^\{x)Hpy^ + X\u{x)Gpp - XupGpx) (34) 

and 

H'pux = Hp^x+dpXv\{.x)+i^{x)Gpx+'n''P{xKv{.x)Hppx-XKv{.x)HpXp-XvK{.x)Hppx+XuK{.x)HpXp). 

(35) 

D. Equations for evolution of matter 

The continuity equations are standard to Yang-Mills theory. By Noether's theorem, for 
any matter action, Sm-, e.g., Dirac's, we have the relation between the potential Gpy and 
the stress-energy-tensor T^'^ and the potential Hp^x and the spin-density S^'^^, 

-■' - wtr '""^ 
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and the continuity equations are, 

^%. = 0, (38) 

D^S^^x = (39) 

[20], where for a generic source J^" and potential A^", D^J,, = d^J^ — i{A^,Jy] is the 
covariant derivative. Hence the full continuity equations are, 

T^xH,^, + T%H^^^)=0, (40) 

5^5M«/3 + G'^^T,^ - G'^f.T^^ + n-^ - n^" = (41) 

where 

na/3 = Tj'^^ {H^aaS^pP — H^^pS^fip — 

H^aaSppi3 + H'^aaS^l3p) ■ (42) 

There are nonlinearities in the field equations 29-30 and the continuity equations because 
gravitons have spin, polarization and momentum and couple to themselves. 

When H^ijx = 0, the field equations without an angular momentum source revert to the 
Abelian equation, 

d'^id^G^x - d,Gf.x) = SttT,,, (43) 

which, under the harmonic gauge condition, is, 

DGp, = SttT^,. (44) 

From this equation, the vortex model, and the geodesic equation derived in the next section, 
I will show that the theory agrees with all observations of N-body motion. 

E. Equation of test particle motion 

First write the SO (4,1) Yang-Mills matrix potential in component form: A^ = GpyV^ + 
H^j_u\M'^^. The test particle has no internal structure and is represented simply as a param- 
eterized path x^{t). Given the basic structure of the theory from Section JIB, the potential, 
G^j,, takes the place of the metric normally found in Lagrangians for test particle motion. 
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Indexes are raised and lowered with the Minkowski metric rjfj^^. The Lagrangian for motion 

is: 

dx^^ dx" dx^ dx^ dx^ , , 

^"^ da da ^^ da da da 

This describes both the curvature and twisting of the path as it moves through spacetime 

along its shortest path. (One can derive this equation by looking at the discrete lattice in 
Section JIB and describing a path as a sequence of connected edges in the lattice.) 

No observations of gravitational phenomena, save perhaps cosmological observations, are 
sensitive enough to measure torsion. Therefore, we concern ourselves in this section only with 
geodesic motion. Let H^^x = 0. (A similar assumption is made in applying Einstein-Cartan 
theory.) Let the Lagrangian for particle motion be, 

^ dx^ dx" ,,„. 

The equations of motion are given by the Euler-Lagrange equation, 

A. ( sc \ _ ar ^ ^ 

da \d{dx^/da) J dx^, 
The equation of motion for the particle (the "geodesic" ) is, 

^ + ^{G-')''' [d,G,. + d,G,, - 9,GJ ^^ = 0, (48) 

which is the geodesic equation. Note: there is a distinction between raising indexes and 
inverting G^^. Its inverse is {G-^Y'' such that {G-^^Gc^u = K^ not G^"" = rj^^'^'q^^Gap. 

F. Strong Equivalence Principle 

Equivalence principles are defining features of gravitational theories. Newton's theory 
included the Weak Equivalence Principle (WEP) where an object's weight was assumed to 
be proportional to its mass. Einstein then developed the equivalence principle named for 
him, the Einstein Equivalence Principle (EEP), which states that local experiments in free 
fall are independent of position and velocity. When strongly gravitating bodies are taken 
into account, the validity of the equivalence principles come into doubt however as internal 
structure may cause violations of WEP such as the (never observed) Nordtvedt effect [1 1] . 

The Strong Equivalence Principle (SEP) is a natural extension of early equivalence prin- 
ciples. It has three conditions: (a) an object's weight is proportional to its mass (WEP) 
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for self-gravitating bodies as well as test bodies, (b) the outcome of any local experiment 
is independent of the velocity of the (freely falling) apparatus, and (c) the outcome of any 
local experiment is independent of where and when in the universe it is performed [11]. 

Although a rigorous proof has not been found, GR is the only known metric theory 
that strictly satisfies SEP. Because experiment has never found a violation of SEP, any new 
gravitational theory must satisfy strict constraints on SEP violations. In the following I show 
that the Yang-Mills theory satisfies SEP strictly (1) because it is not a metric theory and 
its artificial geometry is not affected by gauge transformations and (2) because its gauge 
symmetry allows the required transformations that satisfy SEP as described in the next 
paragraph. 

By satisfying all three conditions, the Yang-Mills theory satisfies SEP. The Yang-Mills 
theory satisfies the first condition, (a), by the arguments of Section 20.6 of [7] also found in 
[21]. If the field of a self-gravitating body asymptotically approaches "flatness" {G^i,{R) = 0) 
at some distance, R, considered to be the boundary of the local system, this is sufficient 
to guarantee that a body's self-gravitation and other internal structure does not affect its 
motion. Because it is always possible to find a gauge that eliminates the gravitational field 
at the boundary between a local, spherically symmetric, compact system and the external 
environment, a spherically symmetric self-gravitating body, even a neutron star or black 
hole, can be regarded as a point particle. In bimetric theories, which are diffeomorphism- 
covariant, changing the gauge to eliminate the field at the boundary changes the Minkowski 
metric, so a coordinate system cannot be found that satisfies this requirement. With the 
Yang-Mills theory, however, which also incorporates the Minkowski metric, the gauge can 
be changed without changing the coordinate system. The Minkowski metric field 1]^,, that 
gives distances in the artificial coordinate system does not change when the gauge of the 
gravitational field G^^, changes. Therefore, the Yang-Mills theory satisfies condition (a). 

The Yang-Mills theory also satisfies the second condition, (b), by local Lorentz covari- 
ance (a subset of local de Sitter covariance) which means that experimental outcomes are 
independent (by gauge covariance) of Lorentz boosts and rotations and satisfies (c) by hav- 
ing no preferred location/time in the field equations. Unlike Rosen's bimetric theory which 
has a local gravitational constant that depends on the field, the Yang-Mills theory has a 
non-location specific gravitational constant, ag. Thus, all three conditions are met, and the 
Yang-Mills theory satisfies SEP. 
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G. Quantization 

The quantization of the theory is given by the generating functional: 

Z[T, S]= f DGDH exp -'- f d^x E^.xE^"^ + F^^^^F^"''^ + la^ f d^xG^.T^^ + H^.^S^""^ 

(49) 

Because it is a massless gauge boson Yang-Mills theory, the theory given in this paper 
has a finite polynomial, positive definite action (20). Since generic Yang-Mills theories on 
semi-simple groups are renormalizable [5], this one is as well. This constitutes mathematical 
proof that the S0(4,l) Yang-Mills theory in this paper is renormalizable. 

The standard quantization scheme is appropriate with Fadeev-Popov ghost fields used to 
derive Feynman rules including the graviton propagator. 



p^ + ie 



(50) 



where \, p are group indexes [5]. 

III. EXPERIMENTS AND OBSERVATIONS 



Numerous observations, starting with light bending in 1919, have been made to attempt to 
confirm predictions of general relativity. These include gravitational time dilation, redshift, 
and light bending which have been measured within the solar system as well as perihelion 
precession of the planets and binary pulsar precession [6] . Other observations do not agree 
with the original theory and have required modifications to the Einstein field equations. 
In order to show that the theory agrees with as many observations as possible, in the 
following, I derive (1) a static, spherically symmetric solution to the field equations, (2) the 
IPN equations of motion for a binary system, (3) the radiation reaction for binary pulsar 
inspiral, and (4) predict the expansion of the universe via a homogeneous, isotropic model. 

A. Spherically symmetric solution 

The static, spherically symmetric Schwarzchild solution is one of the most important 
solutions to the Einstein field equations. The solution to the de Sitter abehan field equations 
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44 is identical up to linear order. Let spacetime be covered by spherical coordinates (t, r, 0, 6) 
where 6 is colatitude and let the metric be the flat spacetime metric [6], 

ds^ = -de + dr"" + r'^de'' + r^ siYi^{e)d^''. (51) 

The only allowable non-zero components of the potentials are Goo and G^r- Choosing a 
harmonic gauge we have, 

V'Goo = SttToo, (52) 

V'^Grr = SllTrr. (53) 

Placing the source at the origin, in spherical coordinates, its tensor has components Too = 
Trr = p such that its density is p = M6{r). Recalling the flat space boundary conditions 
that G^i, — )■ rj^j, as r — )■ oo, the solutions to the field equations are: 

Goo = - (l - V") ' ^- = ^ + V"' ^'' ^ ""'' ^'^'^ ^ ""' '''''^^^- ^^^^ 

This is the linearized Schwarzchild solution to the Einstein equations [6] but, for this Yang- 
Mills field theory, it is exact for all spherically symmetric bodies, including black holes. 

B. Particle trajectories 

Particle trajectories in a static, spherically symmetric potential field are identical to those 
of general relativity. The only difference is in the radial potential. The Schwarzchild radial 
potential grr = (1 — 2M/r)'^ differs from the radial potential derived in Section III A to 
quadratic order, 

Grr = l + 2M/r ^ (1 - 2M/r)-^ + OHM/rf). (55) 

No tests within the Solar System, whether with respect to the Sun, Earth, or another 
body (the most precise tests have been done near the Earth) have achieved better than linear 
order in Schwarzschild coordinates [7] [6] because of the weak fields involved. If M = Gm/c^ 
is the Schwarzchild radius and m is the mass in other units (e.g., kilograms) at the surface 
of the Sun, 2MQ/r ^ 4.25 x 10"^ and at the surface of the Earth 2M^/r ^ 1.4 x 10~^ 
Both are much smaller than necessary to measure higher order effects where imperfections 
in density (e.g., mountains) would make measurements difficult to verify. 
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Predictions of gravitational time dilation and redsliift, meanwhile, are mathematically 
identical in both theories (to any order) given a static, spherically symmetric field with 
redshift proportional to 1 + IMjr. Recent high precision measurements of redshift using 
cesium atoms in a laboratory agree with predictions of both theories [22]. Therefore, no 
tests done within the Solar system to date disprove either theory. For strong field tests, we 
look outside the solar system to binary pulsars and, further on, cosmology. 



C. Post-Newtonian Equations of Motion 

The binary pulsar system B1913+16 discovered in 1974 provided one of the first tests 
of strong field general relativity [21]. Because of the regular radio pulses of the star which 
is orbiting a relatively inert body such as a dead pulsar, henceforth called the companion 
body, allowed extremely precise measurements of the pulsar's orbit, the small deviations of 
the orbit from Kepler's laws can be measured, including orbital damping caused by gravi- 
tational radiation, which is a higher order effect than orbital precession [21]. Observation 
of this system over the past thirty-five years have only increased the accuracy of the mea- 
surements. Recently, the double binary pulsar system PSR J0737-3039A/B with two active 
pulsars orbiting each other has been observed providing even greater confirmation of GR's 
predictions [8] . The following demonstrates that the field equations of Sec. II C also predict 
the observations of these pulsar systems. 

Although it is not a metric theory, its close resemblance to a metric allows the G-field to be 
expressed in a form similar to PPN form. Therefore, in this section, I apply this framework 
to show that the many-body equations of motion agree with those of general relativity up 
to first post-Newtonian corrections. The parameterized post-Newtonian (PPN) equations 
of motion of gravity are the primary vehicle by which not only different theories of gravity 
(including GR) are compared but provide the equations for predicting the motion of two 
or more gravitating bodies. Tests of these equations include the measurement of planetary 
motion within the Solar system for weak fields and binary pulsar precession and gravitational 
radiation for stronger fields [7] [11]. 

For a binary pulsar system the relevant metric in the post-Newtonian coordinate system 
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for general relativity is ([11], 11.52) 

(700 = -l + 2 5^m„/|f-x„(t)|+0(e^), (56) 



a=l,2 
,3\ 



9oj = 0(6^), (57) 

gij = S,ji 1 + 2 Y, mj\x-xa{t)\] + O(e^), (58) 

\ a=l,2 / 

with mi and 1712 being the masses of the bodies and xi and X2 their positions in the appro- 
priate harmonic coordinate system. From the relaxed field equations (44), 



Goo = -1 + 2 V ^^ + 2^ + 0{e'), (59) 



a=l,2 



where 



-r / ^x V^ ^a'i'^ V^ "^a "^h /^^x 

^t,f = V 1^ "^1 + V 1^ °-.||^ ^| , 60 

a=l,2 a,b=l,2,aj^b 

G»i = -2E]|^ + 0(a. (81) 

a=l,2 ' ' 



and 



G,, = 6,,{1 + 2 V — ^) + O(e^). (62) 

a=l,2 1-^ -^"1 

The velocities of the bodies are v\ = (fii,fi2,'yi3) and V2 = ('y2i,'y22,'y23) and Va = \\va\\- 
(See Appendix B for a derivation of these potentials.) 

Dropping the extra terms from G^^, for this system that exceed the required order, the 
IPN metric matches G^^. Since the theory satisfies SEP, for any set of compact nearly- 
spherical bodies where tidal forces may be neglected, the geodesic equation 48 predicts 
orbital motion. The post-Keplerian parameters of periastron advance, (u), time delay, 7', 
and Shapiro delay parameters, r and s, are also, consequently, the same in both theories. 

The final post-Keplerian parameter, orbital speed-up caused by gravitational radiation, 
Pb, derives from the quadrupole formula from the Yang-Mills field equations discussed in 
the following section. 

1. Orbital speed-up of a binary pulsar system 

Gravitational radiation was first addressed by Einstein shortly after the publication of 
general relativity, where he demonstrated that the radiation is primarily quadrupolar in 
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contrast to the primarily dipolar radiation from electromagnetic sources [23]. The energy 
loss of a binary system of compact stars was first demonstrated in a paper by Peters and 
Matthews [24] , who derived the energy loss of binary stars in Keplerian orbit leading to the 
formula for the orbital speed-up, 

M +— e^ + — e*^ j mpmc(mp + mc)"^/^, (63) 

with vrip and vric the masses of the pulsar and companion body respectively and e the orbital 
eccentricity [25]. Although this equation was developed for Keplerian orbits, it has been 
applied to post-Keplerian orbits since the only requirement is that the orbit be elliptical, 
and, although much higher order derivations have been made [11], only this equation has 
been tested. In this section, I demonstrate an identical increase in orbital speed for the 
Yang-Mills theory via the same techniques but omit a general discussion of gravitational 
waves. A detailed discussion of multipole expansions for gravitational radiation in general 
relativity can be found in [26]. 

Orbital speed-up is a function of energy loss, i.e. radiation-reaction, which to lowest 
multipole order in General Relativity is ([7], 36.31), 

dE _ _l ,d Qjk d Qjk , , , 

'dt~~5^^U^ dt^^' ^ ^ 

where Q^^ is the reduced quadrupole moment (the trace free part of the second moment of 
the mass distribution) such that Q^j, = g^j, — \5^uq, 

„.^/d3.p..,..., (65) 

and q = q^^. The Peters-Matthews formula for energy loss of a binary system ([11], 10.80), 

dE I ll^Wb^ 8 / n ox \ , , 

IF = -(Vl5(i^" -"'■■))• ('''" 

where m = mi -|- m^-, fi = 17117712/171, v = \vi — V2\, and r = \xi — X2\, can be found from 
the linearized vacuum equations of general relativity, Dh'^'^ = 0, [24] [7] by the quadrupole 

relation, 

2 
hjkit, x) = -qjkit - r), (67) 

for the spatial part of the radiation field ([7], 36.50) where r is the distance from the source. 
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To show that both GR and the Yang-Mills theory have the same orbital speed-up, we 
compare the equations for the two theories. Peters and Matthews derives the energy loss 
using the linearized Einstein vacuum equations: 

Qfjiu = Vf^u + hf,^, (68) 

such that 



□V = 0, (69) 



where h^^ = /i^^ - \r]^^hx^- 



2 

The linearized Yang-Mills equation for an unpolarized, spinless source, 



UG^, = 8nT^,, (70) 



is the same as the linearized equation for GR. 

Let G^i, — rj^y = h' . A plane wave solution to the wave equation. 



□C = 0, (71) 



IS 



h'^j^ = ae^u cos{ut — k ■ x), (72) 

where a is the amplitude, e^,^ is a symmetric, traceless, transverse, and unitary polarization 
tensor (as defined in [24]), and u is the frequency. 

In the following, I derive the quadrupole expansion for plane wave solutions to the lin- 
earized Yang-Mills (YM) theory. An important difference between electromagnetism and 
gravity is that, whereas the primary electromagnetic radiation for charged bodies in motion 
is dipole, the primary gravitational radiation for massive bodies in motion is quadrupole. 
This derivation follows the outline in [11], Chapter 10. 

The multipole expansion of the plane wave solution (72) is, 

oo ~ 

h'^^ = Ar-^ J]](l/m!)(9/9t)™ / d^x%^{t - r,x'){n ■ x')"'. (73) 

m=0 "^ 

to IPN order, where n = x/r. Because of the gauge condition, d^G^^ = 0, and the retarded 
potential, we can use the following relations, 

doh'ok = nWoh'jf,, doh'oo = n^n'^doh'^j,, (74) 
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and only need to determine the /iL components [11]. Because, to order, d^T^^ = by 
conservation and that the source is symmetric, a useful relation is, 

{d^/de) I d-^ xT^'^x^x'' = 2 fd'^x T^^ (75) 



Therefore, 73 becomes. 



h[j = 2r ^{d^/dt^) I / d'^xT^'^{t — r,x)xiXj j + higher order. 



(76) 



This means that monopole and dipole moments of T*-^ may be expressed as time derivatives 
of quadrupole moments of T°°. To post-Newtonian order T^" = p, the mass density, and we 
now have the quadrupole approximation. 

For the binary system the integral simplifies to, 

Kj = lQAt-r)- (77) 

Since from 67 we have h'^, = hij and because hij solves the vacuum equation we also have 
hij = hij. The speed-up formula (63) follows the same procedure as that of general relativity: 
take a Taylor expansion of 77 in powers of r to extract the radiation-reaction potential and 
applying the geodesic equation 48, determine the reaction acceleration of the bodies ([7], 
Sec. 36.11). Then arrive at 63. 

In this section, I showed that the inspiral of orbiting binary pulsars is the same for both 
general relativity and the Yang-Mills theory. Because the most accurate measurements to 
date match 63, they agree with both theories. 

D. Cosmology 

If this paper had been written as recently as 10-15 years ago, cosmology could have 
been addressed more fully within the scope of this paper. In recent years, however, the 
amount of data available to match to cosmological models has exploded thanks to projects 
such as WMAP [4], observations of Type la supernovas, and baryon acoustic oscillation 
measurements. In this section, I briefly address the Robertson- Walker model, the starting 
point of the two most prominent numerical approaches, both centered around the Friedmann 
equations: perturbations of the Einstein field equations and N-body cosmology which models 
the universe as an N-body Newtonian system. 
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The Friedmann equations can be derived from Newton's laws, suggesting that at slow 
speeds and very weak fields they are valid for the universe as a whole — certainly the present 
universe is essentially Newtonian at large scales where inhomogeneities are small. In the 
early universe, however, relativistic speeds and strong gravitational fields predominated and 
the Newtonian model breaks down for times "close" to the Big Bang. 

1. Isotropic, Homogeneous Universe 

Observations of the cosmic microwave background and statistical counts of the distri- 
bution of galaxies imply that the universe is statistically isotropic [6]. By the Copernican 
principle that we do not occupy a privileged location in the universe, the universe must also 
be homogeneous. An isotropic, homogeneous universe has a very simple dynamical descrip- 
tion. For a perfect fiuid universe, the stress-energy-momentum tensor in special relativity 
has the standard form, 

T^- = pu^'u'' + piu^'u'' + g^"") , (78) 

where p is mass density and p is pressure [6]. The fundamental Friedmann equations for 
general relativity are, 

3a 



,2 



a' + k = :^, (79) 



-87r(p + 3p), (80) 

where k is the curvature constant [10] and 

5-00 = -1, grr = ^ _ , 2 , gee = a'^{r)r'^, g<t><t> = a^{r)r^ sin^ 0, (81) 



is the metric. Recent observations that the universe's density, the sum of baryon, dark 
matter, and dark energy densities, is near or at the critical value, Qb + ^c + ^A — ^crit, and 
the universe is spatially fiat, k ^ [4] [9]. Therefore, for the rest of this section. A; = 0. From 
the conservation of energy, 

|^ + -(P + P) = 0, (82) 

OT a 

is a fundamental equation for the evolution of matter and energy. 

Turning to the Yang-Mills theory, I derive the equations of motion for the universe. 
Modern cosmology relies heavily on perturbation methods, and the Robertson- Walker as- 
sumptions of isotropy and homogeneity are starting points for this approach because they 
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generate the unperturbed solutions to the field equations. Therefore, I again make these 
assumptions. 

The isotropic, homogeneous universe under the General Relativity model does not dis- 
tinguish between matter moving apart within space and space itself moving apart. Hence, 
redshifts are interpreted as being caused by matter moving apart as space drags it along in 
its expansion. It is certainly possible, however, to distinguish between the two types of red- 
shift, gravitational and doppler, in ordinary, non-cosmological circumstances such as within 
the Solar System. Making such a distinction is important when attempting to measure (1) 
the age of the universe, which is based on redshift measurements, and (2) whether matter 
is accelerating away or experiencing increasing gravitational redshift and moving away at 
fixed velocity. These measurements, in turn, affect how observations such as abundances of 
light elements are interpreted [16]. 

In this section, I show how the theory supports a linear coasting or Milne cosmology with 
some modifications: In the Poincare approximation, gravitational redshift is linear. Thus, 
in that approximation all observed redshift is caused by matter fiying apart. In the de Sitter 
theory, however, the infiuence of the torsion potential, H , causes gravitational redshift to 
increase causing redshift to accelerate with time and giving the impression of an accelerating 
expansion. 

As discussed in section II B, the coordinate system is a set of degrees of freedom, that, 
while having no physical meaning, may be used to simplify equations. In this cosmology, I 
will use comoving coordinates again, g^^ = —l,gij = a^{t)6ij, to expand with matter moving 
within space and define the expansion of space itself with the following. Goo = ~1 smd 
Gij = b{t)a^{t)6ij. The gauge is also comoving. As we will see, the torsion potential is 
non-zero, -f^^i^A 7^ 0, because where two G gravitons couple they produce H gravitons. This 
is a consequence of the de Sitter Lie algebra [Vfj_, Vy\ ^ and disappears in the Poincare 
approximation where [V^, Vy\ = (see B). Let i/joj = —Huq = c{t)a'^{t) and Huj = —Hiji = 
d{t)a^{t) be the isotropic, homogeneous torsion potential. 

In the comoving coordinates, the conserved current is Too = P and Tij = pa^5ij. The spin 
density is zero by isotropy, S^j^yx = 0. Evaluating the field equations 29, we have, 

3a(ab + ba — a + ca) , . 

— i = Snp, 83 
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and the next three (all the same) are, 

—3aba — ha^ — 'dab — a% + ad ~ ha^ — ca^ — Idac 



a2 



iTiV- (84) 



From the field equations 30, we have, 

ca^ + 3cda — 3d^c + dac — dab + ad — ba^ — ca^ = (85) 

and 

da^ — a^d + add — a^b'^d = 0. (86) 

By the vortex model all matter has p = p/3 which simplifies 40 to 

P = -4-p (87) 

a 

The system of three differential equations (83,84, and 87) govern the gravitational behavior 
of matter in this universe. The boundary condition are given in terms of present day, to, 
and origin, t = 0, parameters: a(to) = Oq, a(0) = 0, 6(0) = 0, c(0) = 0, and c(to) = cq. The 
assumption of a(0) = is a classical assumption. Quantum effects take over at t ^ 1 and 
may predict a different initial size of the universe. Certainly, the matter in the universe could 
never be localized to a point of zero radius as this would violate the uncertainty principle. 
The equations 83,84, and 87 have the non-trivial solution, 

a{t) = %, (88) 

to 

Pit) = Po4' (89) 

^(^) = Y^ [^^' - 3a-^ + (12 + Snpott)^ - AOnpo4t - 3271 po4] , (90) 

c(t) = ^[-3Ct + 87rpo4], (91) 

where C = -^ (— 12co + SnpotQ). The equation for d{t) is not expressible in closed form but 
is the solution to a second order ODE. 

The rate of expansion of the universe, v, is related to the Hubble constant by, 

y^^ = HR, (92) 

dr 

where R is the distance between two isotropic observers [6], and 

H =- = lit, (93) 

a 
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where currently accepted value of the Hubble constant is 70.1 ± 1.3{km/s)/Mpc [4]. This 
is the apparent Hubble constant, however. Let /3(t) = ^/b{t). Then the apparent Hubble 
constant is: 

H=^^ = ^^+H. (94) 

pa p 

This implies that the universe's age should be given by, 

^0 = -^. (95) 

If 4 > the universe is older than it appears. If 4 < it is younger than it appears. Given 
that the redshift is accelerating, the former is true, and the universe is older than it appears. 
This makes sense because gravitational redshift on top of Doppler redshift would make the 
expansion of the universe appear faster than it is implying that its origin was sooner in the 
past than the Doppler redshift alone would imply. 

The scale factor, a(t), conforms to a universe with a Minkowski metric sometimes called 
the Milne universe or simply the linear model [27]. The Milne cosmological model has 
recently seen a resurgence as an alternative to the standard ACDM model and has some 
compelling features: It gives the same age for the universe as the ACDM model [28]. Nu- 
cleosynthesis has been fit to observations for the linear model in [29] with some problems 
that I discuss below [16]. Predictions for Type la supernovae are very close to those of the 
standard model. The Milne universe has angular distances in the CMB of about 1.1 degrees, 
also close to the observed value [15] [30]. 

A variation known as the "symmetric" Milne model makes use of equal amounts of matter 
and antimatter present in the universe to account for zero energy density, assuming that 
antimatter has negative mass [15]. This assumption is not necessary here. Indeed, antimatter 
in the Yang-Mills theory has positive mass and energy (Section II A). 

Although the rate of expansion is linear, the gravitational potential, 6(t), shows an ac- 
celeration in apparent expansion rate, related to y^G^j^/S = ^Jh{t)a?{t) ~ 0{t'^), meaning 
that the redshift of distant galaxies should be accelerating linearly with time even while the 
galaxies are actually moving away at a constant rate given by a{t). At the present time, the 
acceleration is not well understood enough to determine whether it is linear, exponential 
(as in models with a cosmological constant), or some other function of time. Hence, should 
better measurements become available that show the rate of acceleration, they may test the 
theory's predictions. 
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2. The Poincare approximation 

When the Poincare approximation is made such that we change the equations to conform 
to the Lie algebra with [V, V] ~ (see B), the expansion no longer accelerates: 

a(t) = ttot/to, (96) 

p{t)=po{^)\ (97) 

m = 1 - TTTT^ + r^TV,^ 98 

Siayf^jV {ao/to)t 

where K = (6o + 87rpo/(3ao/^o) ~ l)'^o and bo = b{to). Note how b(t) — )> 1 as t — )> oo 
indicating that the Poincare universe has no accelerating expansion. The explanation is 
that when [V^, V] = iAf^" (see A), it means that momentum, related to the V^ generators, 
produces spin. Because [M^^,Pp] = iirj^pVi, — rj^pV^) spin gravitons couple to momentum 
gravitons to create more momentum gravitons. Since the potential for momentum gravitons 
is G^y this relationship causes the redshift we see. The result is an accelerating expansion 
of gravitational redshift. 

Another physical phenomenon has a similar explanation: Thomas Precession. In that 
case, the Lorentz rather than de Sitter algebra works to relate orbital momentum of the 
electron to changes in its spin [14] . Thus, "dark energy" can be explained as the production 
of gravitons by an interplay between spin and momentum. 

IV. POTENTIAL DISAGREEMENTS 

Observations of gravity can be classed into a four basic types: (1) Newtonian gravity 
which applies to very weak fields and slow velocities that pertain to most galaxies, stars, 
and orbits of the outer planets, (2) weak field gravity responsible for classic effects such as 
redshift, light bending, time dilation, perihelion precession (3) strong field gravity responsible 
for binary pulsar precession, orbital damping, and gravitational waves, and (4) cosmology 
where assumptions such as homogeneity and isotropy greatly simplify equations as well as 
perturbative cosmology (not addressed in this paper). Included are tests of the Strong 
Equivalence Principle (SEP). 

Going through each of the four types of tests: 

31 



Newtonian gravity. Both theories subsume Newtonian gravity at very weak field strengths 
and slow velocities, and this is easily demonstrated from the field equations. 

Weak Fields. As shown in section III B, both theories agree with classic effects up to linear 
order in the Schwarzschild radius over the distance. Measurements of redshift, light bend- 
ing, perihelion precession, and time dilation have not been made to quadratic accuracy in 
Schwarzschild coordinates (as opposed to isotropic coordinates where the coordinate change 
introduces an artificial "measurable" quadratic term); therefore, none of these experiments 
contradict the YM theory. 

Strong Fields. As [8] mentions, measurements of binary pulsar precession, while some of 
the most precise measurements of relativistic gravity ever made, are not sufficiently accurate 
to confirm beyond the IPN equations of motion. Orbital damping caused by radiation 
reaction is a higher order effect but also not confirmed beyond the quadrupole approximation, 
and neither GR nor the YM theory predict a dipole moment. Therefore, as shown in section 
IIIC, none of these measurements contradicts either theory nor any other measurements of 
multibody systems. The YM theory's geodesic equation 48 is also the same as in relativity; 
therefore, since SEP is satisfied (Sec. II F) and compact self-gravitating bodies behave as 
if they are test bodies, it agrees with the IPN equations of motion of general relativity. 
Gravitational wave measurements promise to provide higher order estimates which may 
show a violation of one of the theories, but these, as yet, are not available. 

Cosmology. The highest order measurements available are cosmological. In this case, 
several sources including Type la supernovae (SN), baryon acoustic oscillation (BAO), and 
the cosmic microwave background (CMB) provide data that can be used to constrain various 
theories [9]. The model we have presented is essentially a toy model because it assumes no 
fiuctuations or deviations from perfect uniformity and isotropy in the universe. First order 
perturbation of the field equations about the Robertson- Walker solution combined with 
numerical solvers can indicate whether observed fiuctuations in the CMB agree with the 
theory, and this is an important future direction of research for the YM theory. 

Studies of Big Bang Nucleosynthesis (BBN) where the predicted quantities of light el- 
ements such as Helium, Lithium, and Deuterium are compared with predicted quantities 
can also constrain some theories and several studies of Tensor-Scalar theory have been done 
[31] [32] . These studies tend to focus on the effects of the scalar field on nucleosynthesis and 
are not directly applicable to the YM theory which has no scalar field. 
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While the hnear model explains features of the later universe as well as the standard 
ACDM model [15], a single power law scale factor has proved insufficient to explain the early 
universe where the power is sharply constrained to around a ~ t°'^^ [16] unless additional 
modes of production of deuterium are invoked in the later universe. Because the YM theory 
predicts a gravitational redshift on top of the Doppler redshift, however, the age of the 
universe may be different than currently predicted. Since the age is an essential assumption 
of nucleosynthesis predictions, this throws the arguments against the linear model, in the 
context of the YM theory, into doubt. Whether the theory is able to predict abundances of 
light elements is left for future work. 

V. FUTURE TESTS 

General Relativity and the Yang- Mills theory agree for the linearized Schwarzschild metric 
predictions and IPN equations of motion but, with the exception of radiation-back-reaction 
where the same degree of inspiral is predicted for binary pulsar systems, they fundamentally 
diverge at higher order motion. The 2PN order should be measurable for binary pulsars 
within the next few decades as additional data is collected about known systems (the more 
data collected the smaller the error bound). Gravitational wave detectors such as LIGO and 
VIRGO may also be capable of detecting much higher order effects (well into the quadratic 
or even cubic realm) from inspiraling black holes and other closely spaced, highly massive 
objects. Because general relativity and the Yang-Mills theory no longer agree at this level, 
these observations will determine whether gravity is truly diffeomorphic or a non-Abelian 
gauge theory like the other known forces. Another important test is a quantitative mea- 
surement of the values of h^ and Cq which may allow a quantitative prediction of the rate of 
acceleration of the universe, an important prediction of the theory. 

Quantum predictions are significant because the theory, in order to be renormalizable, 
requires mass, energy, momentum, and angular momentum to be quantized so that the 
coupling constant has zero mass dimension. Without that assumption, the coupling constant 
would have negative mass dimension of -2 and the theory would not be renormalizable [5] . 
The size of the quantum must be exceptionally small, but an experiment may be devised 
to detect it similar to how Millikan's oil drop experiment detected the indivisibility of the 
electron charge. 
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VI. RELATED WORK 



While this paper has focused on the problems with the prevailing theory of gravity, general 
relativity, several other theories have been proposed to solve quantum gravity and explain 
macroscopic observations. The most important empirical theory is the Lambda-CDM model 
which was used in the WMAP survey [4] and combines the Einstein field equations with a 
cosmological constant and cold dark matter. The field equations 29-30 do away with the 
cosmological constant and much of the complexity of the curved spacetime model, but the 
YM theory does not challenge the prevailing theory of dark matter. 

The most prominent theory of quantum gravity is string theory and its derivative M- 
theories which, rather than being a simple theory of gravity, is an extensive modification 
of physical law positing that all matter is composed of strings, attempting to unify all 
forces [5] [33]. Although unification is not its purpose, in representing gravity identically 
to the other three forces, the theory may possibly be unified with them at a high enough 
symmetry group without requiring any significant additional physical assumptions such as 
strings or additional spacetime dimensions. 

The next most significant model of quantum gravity is loop quantum gravity [33]. Loop 
quantum gravity has a different approach from that given in this paper, further marrying 
gravity to curved spacetime. In making the assumption that gravity has only diffeomorphism 
covariance and is dominated by a metric geometry, it retains the complexity of general 
relativity and adds to it by introducing spin foams, i.e. discretizations of spacetime, in 
order to avoid blow-ups of the quantum variables. It also makes no predictions that are 
currently testable and, furthermore, has not been shown to agree with classical general 
relativity in its entirety. 

The main Yang-Mills approach to gravity is conformal gravity [12], and is not directly 
relevant to the theory in this paper which makes significantly different physical and mathe- 
matical assumptions. Conformal gravity is a four derivative Lagrangian theory and predicts 
significantly different large scale behavior than either general relativity or Newtonian gravity. 
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VII. CONCLUSION 

The assumption of curved spacetime, while elegant, is not a necessary component given 
current experimental evidence, and, therefore, we are motivated to consider alternatives to it 
in a theory of gravity. By assuming a de Sitter symmetry group, this theory not only agrees 
with the observations confirming general relativity's major predictions but also includes a 
renormalizable quantum theory of gravity and a definition of mass as relativistic energy. 
It does this by establishing gravity on the theoretical foundation of the Standard Model. 
Indeed, the Standard Model SU{3) x SU{2) x f/(l) is conjectured to be part of a unified, 
covariant theory in SU{5) at some energy since SU{5) is the smallest group that can contain 
all three. With the addition of the de Sitter group, the Standard Model could be updated to 
be 50(4, 1) X SU{3) x SU{2) x f/(l). The smallest group to contain this would be SU{10) 
with 99 degrees of freedom. Another possibility would be 5*0(15) which has 105 degrees of 
freedom. Either of these might represent a unified field theory at high enough energy. 

The theory, in addition, solves at least two problems with General Relativity: 

1. Quantum renormalizability. General Relativity has a non- renormalizable action while 
the YM theory has a renormalizable action provided mass and momentum are quan- 
tized. Although such quantization has not been measured, the weakness of gravity 
suggests that it is extremely small, much smaller than the quantization of charge. 

2. Vacuum energy. Although vacuum energy is the predominant explanation for the 
source of dark energy, the vacuum energy has a density approaching h [5] . This density 
is the density the universe would have if Planck mass black holes were set adjacent 
to one another, event horizon to event horizon at every point in space. Clearly, this 
is not so. It has been suggested that something cancels the vacuum energy out, but 
that has not yet been found. It is reasonable to believe that the coupling of gravity 
to the vacuum energy, far from being an explanation for dark energy, is, along with 
the accelerating expansion, an unsolved anomaly of the mainstream theory. By not 
coupling to the vacuum energy in its quantum formalism, the YM theory does not 
have the same problem. 

The most important prediction, and the reason why this paper is more than a theoretical 
exercise, is the prediction of the accelerating expansion of the universe from the de Sitter 
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group. Although only qualitative here, it has the potential to be a quantitative prediction 
with additional data to constrain the current conditions of the universe within the context 
of this theory, particularly values of 60 and Cq. Further confirmations may be obtained from 
detection of gravitational waves, observations of supermassive objects such as black holes, 
and additional cosmological measurements. Additional numerical studies are required to 
determine how Big Bang Nucleosynthesis observations fit the theory. 
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Appendix A: The de Sitter Group Lie Algebra 

In addition to rotations and boosts, all physical theories are translation invariant. This 
confirms the assumption that there is no "special" place in the universe. By Noether's 
theorem, translation invariance causes momentum to be conserved (including energy). A 
translation by an amount v^ is achieved by a 5x5 matrix, T = v'^Vfj, such that V^ = {Vt, Vj) 
are four generators. Given a 4-vector u^iiw = {u^, 1), then w' = Tw = {u^+Vfj,, 1). It is well 
known, however, that prior to the introduction of Lorentz covariance, mechanical theories 
such as Newtonian gravity were Galilean invariant, i?^®S0(2) rather than Lorentz invariant, 
S0(3,l). A reasonable extension of the Poincare group, then, is from i?*^^'^)(8>SO(3,l) to 
S0(4,l) assuming that the radius of curvature in the fifth dimension is large so that it 
appears to be Poincare invariant for small rotations. 

In the de Sitter group, X4 is a fourth spatial dimension and has ordinary rotations with 
respect to the other three spatial dimensions and boosts with respect to time. The anti-de 
Sitter group has x^ as a second time dimension. Let V^ be rotations/boosts in the x^-X/^ 
plane. 

The Poincare group has the covariant Lie algebra, 

K.vj\ = (Al) 

[M^„ Pp] = i{7]^pV, - r],pVp) (A2) 

[Mp,, Mp„] = i{r]ppM,^ - r]^Mup - VupMp^ + r],Mpp), (A3) 

for rjp^iy the Minkowski metric where U = exp[— |ci;^^M^'^] is a Lorentz operation and U = 
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exp[iafj,V^] is the translation operation [34]. The only alteration that the de Sitter group 
makes is that the "translation" operators now generate a Lorentz rotation or boost, 

[V^,V,] = tM^,, (A4) 

in which case 

is the de Sitter operation. 



Appendix B: Derivation of the IPN potential 

The potential at a spacetime point (t, x) in the given inertial frame can be solved by the 
integral equations [21], /i^^ = (G^,, - r]^^)/2: 

K.{t,x) = /dV ^^-^^-J^-f^''^^^ (Bl) 

Jc \x-x'\ 

where C is the past-light cone and bodies have small polarized acceleration and angular 
momentum is small. 

The integral equation may be solved by iteration in which a trial potential h^^ = is 
inserted into the integrand and the integral solved to arrive at a new potential h' which is 
then reinserted to achieve a potential h" and so on until the required accuracy is achieved. 
In order to derive the IPN equations of motion general relativity requires two iterations but 
here we only require one, meaning that the IPN equations for this theory are linear. 

Let the baryon density be po; and the Newtonian potential is given by, 

f/=/dV?^. (B2) 

J \x — x'\ 

Let f/ ~ e^ and v/c ~ e where v is the average velocity and c is the speed of light in 
vacuum. The first post-Newtonian corrections require that /iqo be given up through order 
e^, hoj through order e^ and hij through e^, correcting the Newtonian order, 

/loo = U, hoj = 0, hij = 0. (B3) 

The G-field derives from the equations of Sec. TIC above and the stress-energy- 
momentum tensor for baryon dust. Each baryon particle or compact spherical mass (such 
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as a neutron star or black hole) p has a stress-energy-momentum tensor in its rest frame in 
a spherical coordinate system with the particle at the origin, 

Too = 'mp5^{x), Trr = mp5^{x). (B4) 

The latter is related to the potential hrr which has a transformation, 

which is 

[Xi — X^)[Xj — Xj) /-r,n\ 

hij = rup — — —:- (B6) 

\x - x']-^ 

for a particle at x'. To obtain the potential for all the dust, sum all the potentials for all 
the baryons by superposition of the potentials, 

hij = / d-'x — — — r '— + 0(e ). (B7) 

J \x — x'\-^ 

By a non-Abelian gauge transformation [7], 



hij = Gij + diOj, (B8) 



X= (rx'po{t,x')\x-x'\, (B9) 



where 9fj, = d^x and 



we have 

hij = SijU + 0{e^). (BIO) 

The next two parts of the G-field, Hqj and /loo, require particle velocities and time dilation 
to be included. The current tensors, T^y{p) of each baryon p, are in the rest frame, but, 
because the bodies are moving with respect to the observer at infinity and subject to the 
gravitational fields of one another, the tensors must be (1) boosted from comoving frame 
x^ to the observer's frame, x^, and (2) subject to a non-Abelian gauge transformation from 
the zero gravitational field (free fall) in the rest frame of the body to the gravitational field 
in the rest frame of the observer, 

dx^ dx^ 
Ti,v{v) = ^^^"/^(^) + h^^mp6^{x-Xp), (Bll) 

where A is a Lorentz transformation matrix for a boost density v{t,x) (c = 1)[7] and 
h^i, is the G-field not including the contribution of the mass at that point. These two 
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transformations are equivalent to the two step, boost-and-coordinate-transform method of 
general relativity [7] which reflects the effects of both velocity and gravitational fields on 
measurements of stress-energy-momentum tensors. 

For a baryon p with mass rrip at x' with a rest frame potential hfj_i, and velocity v {c = 1), 



where 

mp{xi -x'i){xj -x'j) 

\x — x\ 
To order this becomes, 

m„ ( \{x - f) • v\{xj - x'j) 

I rf rft I \ "^ 'T* 'T*' 



% = ^^1.^:^13 ■ (B13) 



^0^- = -^3^ -^- + ^^ 1;_^|2 ^ + O(e^). (B14) 



Under the same gauge transformation as above, 

hoj = hoj -do9j, (B15) 

and summing over all baryons, the potential is, 

V = -2\/,+0(e^), (B16) 



where 



The time-time field is. 



where 



and 



1^,= /dV^^^M>#^. (BIT) 

' \x - x'\ 



/ioo = f/ + 2^ + <l> + 0(e^), (B18) 



nt,x)= I a^,' P2MM±E1^ (B19) 

X — x'\ 



^(M)-/dV§i^| l'V-^2,;f -.n^ (B20) 



The gauge transformation, /iqo = Gqo + d^x (which we have already applied to the other 
potentials and hence must carry over to this one), eliminates the last term, and gives 

/ioo = t/ + 2^ + 0(e^). (B21) 

Missing is a term of order f/^ that appears in metric PPN formalism which indicates the 
nonlinearity in the superposition of the potential. In GR, however, this term is an artifact of 
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the choice of isotropic coordinates and not a true nonhnearity such as would be associated 
with massive gravitons. 

In general relativity, the change of coordinates, r = f(l + M/2f)^, switches the 
Schwarzschild solution to isotropic coordinates [7], 

d-s" = - ( \~ ^J'^^' \ de + {l + M/2ff[df^ + f^{de^ + sm^edcf)]. (B22) 

Because the change in coordinates, however, depends on the potential itself, this introduces 
a nonhnearity that previously did not exist in the solution. Like the Schwarzschild solution 
in Schwarzschild coordinates, the Yang-Mills theory IPN solution does not contain any 
nonhnearity because it is not in this nonlinear coordinate system. The exact, Yang-Mills 
spherically symmetric solution given in spherical coordinates, 

/ioQ = h„ = M/r, (B23) 

under the same change of variables, becomes to IPN order, 

hoo = M/f - 2{M/ff + 0{e^), (B24) 

h^^ = h,o = 0{e') (B25) 

and 

hij = 6ijM/f + 0{e^). (B26) 

While one could make a change of the flat spacetime coordinates to bring the YM potential 
in line with the GR metric, there is no purpose in doing so. Since both GR and the Yang- 
Mills theory are generally covariant, changing coordinates does not change predictions. It 
does, however, require that observations be input in the appropriate coordinate system since 
the Schwarzschild r is not the same as the IPN f. Indeed, because distances are defined 
by the metric, g^^, for the YM theory r is the true distance between a point and the origin 
while, in General Relativity, the Schwarzschild r is simply a coordinate. Meanwhile, because 
of the potential dependent metric in "isotropic" coordinates, f is just a coordinate in both 
theories and does not refer to distance directly. 

The remaining equations (which can be derived from the field equations) are the equation 
of motion. The continuity equations (38) to IPN order are the four equations, 

S^T*^" = 0, (B27) 
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each of which is the same as the electromagnetic continuity equation. Linearized general 
relativity has the same problem [6] and care needs to be taken that the nonlinear, gauge 
covariant equations of motion (38) be used instead. For A^-body problems, which are gener- 
ally the ones that gravitation deals with, however, the geodesic equation is sufficient. This 
is because the Yang-Mills theory satisfies the Strong Equivalence Principle (SEP) and inter- 
nal structure does not affect the motion of bodies. Therefore, any spherical, compact body 
moves as a test particle. 
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